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Abstract. Over the past decade, supernovae have emerged as some of the most power- 
ful tools for measuring extragalactic distances. A well developed physical understanding 
of type II supernovae allow them to be used to measure distances independent of the 
extragalactic distance scale. Type la supernovae are empirical tools whose precision and 
intrinsic brightness make them sensitive probes of the cosmological expansion. Both 
types of supernovae are consistent with a Hubble Constant within ~10% of Ho ~ 70 
km s~^ Mpc~^. Two teams have used type la supernovae to trace the expansion of the 
Universe to a look-back time more than 60% of the age of the Universe. These observa- 
tions show an accelerating Universe which is currently best explained by a cosmological 
constant or other form of dark energy with an equation of state near w = p/ p = —1. 
While there are many possible remaining systematic effects, none appears large enough 
to challenge these current results. Future experiments are planned to better charac- 
terize the equation of state of the dark energy leading to the observed acceleration 
by observing hundreds or even thousands of objects. These experiments will need to 
carefully control systematic errors to ensure future conclusions are not dominated by 
effects unrelated to cosmology. 



1 Introduction 

Understanding the global history of the Universe is a fundamental goal of cosmol- 
ogy. One of the conceptually simplest tests in the repertoire of the cosmologist is 
observing how a standard candle dims as a function of redshift. The nearby Uni- 
verse provides the current rate of expansion, and with more distant objects it is 
possible to start seeing the varied effects of cosmic curvature and the Universe's 
expansion history (usually expressed as the rate of acceleration/deceleration). 
Over the past several decades a paradigm for understanding the global properties 
of the Universe has emerged based on General Relativity with the assumption of 
a homogeneous and isotropic Universe. The relevant constants in this model are 
the Hubble constant (or current rate of cosmic expansion), the relative fractions 
of species of matter that contribute to the energy density of the Universe, and 
these species' equation of state. 

Early luminosity distance investigations used the brightest objects available 
for measuring distance ~ bright galaxies [3,39], but these efforts were hampered 
by the impreciseness of the distance indicators and the changing properties of 
the distance indicators as a function of look back time. Although many other 
methods for measuring the global curvature and cosmic deceleration exist (see, 
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e.g., [66]), supernovae (SNe) have emerged as one of the preeminent distance 
methods due to their significant intrinsic brightness (which allows them to be 
observable in the distant Universe), ubiquity (they are visible in both the nearby 
and distant Universe), and their precision (type la SNe provide distances that 
have a precision of approximately 8%). 

2 Supernovae as Distance Indicators 

2.1 Type II Supernovae and the Expanding Photosphere Method 

Massive stars come in a wide variety of luminosities and sizes and would seem- 
ingly not be useful objects for making distance measurements under the standard 
candle assumption. However, from a radiative transfer standpoint these objects 
are relatively simple and can be modeled with sufficient accuracy to measure dis- 
tances to approximately 10%. The expanding photosphere method (EPM), was 
developed by Kirshner and Kwan [44], and implemented on a large number of 
objects by Schmidt et al. [86] after considerable improvement in the theoretical 
understanding of type II SN (SNII) atmospheres [15,16,99]. 
EPM assumes that SNII radiate as dilute blackbodies 



where 9ph is the angular size of the photosphere of the SN, Rph is the radius of 
the photosphere, D is the distance to the SN, Fx is the observed flux density of 
the SN, and B\(T) is the Planck function at a temperature T. Since SNII are 
not perfect blackbodies, we include a correction factor, which is calculated 
from radiate transfer models of SNII. SNe freely expand, and 



where Vph is the observed velocity of material at the position of the photosphere, 
and t is the time elapsed since the time of explosion, to- For most stars, the stellar 
radius ,i?o, at the time of explosion is negligible, and Eqs. (1-2) can be combined 
to yield 



By observing a SNII at several epochs, measuring the flux density and tem- 
perature of the SN (via broad band photometry) and Vph from the minima of 
the weakest lines in the SN spectrum, we can solve simultaneously for the time 
of explosion and distance to the SNII. The key to successfully measuring dis- 
tances via EPM is an accurate calculation of C{T). Requisite calculations were 
performed by Eastman et al. [16] but, unfortunately, no other calculations of 
C(T) have yet been published for typical SNIIP progenitors. 




(1) 



Rph — Vph{t — to) + Rq, 



(2) 
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Hamuy et al. [34] and Leonard et al. [52] have measured the distances to 

SN1999cm, and they have investigated other aspects of EPM. Haniuy et al. [34] 
challenged the prescription of measuring velocities from the minima of weak 
lines and developed a framework of cross correlating spectra with synthesized 
spectra to estimate the velocity of material at the photosphere. This different 
prescription does lead to small systematic differences in estimated velocity using 
weak lines but, provided the modeled spectra are good representations of real 
objects, this method should be more correct. At present, a revision of the EPM 
distance scale using this method of estimating Vph has not been made. 

Leonard et al. [51] have obtained spectropolarimetry of SN1999em at many 
epochs and see polarization intrinsic to the SN which is consistent with the SN 
have asymmetries of 10 — 20%. Asymmetries at this level are found in most SNII 
[101], and may ultimately limit the accuracy EPM can achieve on a single object 
(cr ~ 10%). However, the mean of all SNII distances should remain unbiased. 

Type II SNe have played an important role in measuring the Hubble constant 
independent of the rest of the extragalactic distance scale. In the next decade 
it is quite likely that surveys will begin to turn up significant numbers of these 
objects at z ~ 0.5 and, therefore, the possibility exists that SNII will be able to 
make a contribution to the measurement of cosmological parameters beyond the 
Hubble Constant. Since SNII do not have the precision of the SNIa (next section) 
and are significantly harder to measure, they will not replace the SNIa but will 
remain an independent class of objects which have the potential to confirm the 
interesting results that have emerged from the SNIa studies. 

2.2 Type la Supernovae as Standcirdized Ccindles 

SNIa have been used as extragalactic distance indicators since Kowal [42] first 
published his Hubble diagram (cr = 0.6 mag) for type I SNe. We now recognize 
that the old type I SNe spectroscopic class is comprised of two distinct physical 
entities: SNIb/c which are massive stars that undergo core collapse (or in some 
rare cases might undergo a thermonuclear detonation in their cores) after los- 
ing their hydrogen atmospheres, and SNIa which are most likely thermonuclear 
explosions of white dwarfs. In the mid-1980s it was recognized that studies of 
the type I SN sample had been confused by these similar appearing SNe, which 
were henceforth classified as type lb [59,94,102] and type Ic [36]. By the late 
1980s/early 1990s, a strong case was being made that the vast majority of the 
true type la SNe had strikingly similar light curve shapes [11,46-48], spectral 
time series [6,18,28,62], and absolute magnitudes [47,54]. There were a small 
minority of clearly peculiar type la SNe (e.g., SN1986G [63], SN1991bg [19,49], 
and SN1991T [19,78]), but these could be identified and removed by their un- 
usual spectral features. A 1992 review by Branch and Tammann [7] of a variety 
of studies in the literature concluded that the intrinsic dispersion in B and V 
maximum for type la SNe must be < 0.25 mag, making them "the best standard 
candles known so far." 

In fact, the Branch and Tammann review indicated that the magnitude dis- 
persion was probably even smaller, but the measurement uncertainties in the 
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available datasets were too large to tell. The Calan/Tololo Supernova Search 

(CTSS), a program begun by Hamuy ct al. [31] in 1990, took the field a dra- 
matic step forward by obtaining a crucial set of high quality SN light curves 
and spectra. By targeting a magnitude range that would discover type la SNe 
in the redshift range z = 0.01 — 0.1, the CTSS was able to compare the peak 
magnitudes of SNe whose relative distance could be deduced from their Hubble 
velocities. 

The CTSS observed some 25 fields (out of a total sample of 45 fields) twice 
a month for over three and one half years with photographic plates or film at 
the Cerro Tololo Inter-American Observatory (CTIO) Curtis Schmidt telescope, 
and then organized extensive follow-up photometry campaigns primarily on the 
CTIO 0.9 m telescope, and spectroscopic observation on either the CTIO 4 m 
or 1.5 m telescope. Toward the end of this search, Hamuy et al. [31] pointed 
out the difficulty of this comprehensive project: "Unfortunately, the appearance 
of a SN is not predictable. As a consequence of this we cannot schedule the 
foUowup observations a priori, and we generally have to rely on someone else's 
telescope time. This makes the execution of this project somewhat difficult." 
Despite these challenges, the search was a major success; with the cooperation 
of many visiting CTIO astronomers and CTIO staff, it contributed 30 new type 
la SN light curves to the pool [32] with an almost unprecedented control of 
measurement uncertainties. 

As the CTSS data began to become available, several methods were presented 
that could select for the "most standard" subset of the type la standard candles, 
a subset which remained the dominant majority of the ever-growing sample [8]. 
For example, Vaughan et al. [97] presented a cut on the B-V color at maximum 
that would select what were later called the "Branch Normal" SNIa, with an 
observed dispersion of less than 0.25 mag. 

Phillips [64] found a tight correlation between the rate at which a type la 
SN's luminosity declines and its absolute magnitude, a relation which apparently 
applied not only to the Branch Normal type la SNe, but also to the peculiar type 
la SNe. Phillips plotted the absolute magnitude of the existing set of nearby 
SNIa, which had dense photoelectric or CCD coverage, versus the parameter 
Z\mi5(B), the amount the SN decreased in brightness in the B-band over the 15 
days following maximum light. The sample showed a strong correlation which, 
if removed, dramatically improved the predictive power of SNIa. Hamiiy et al. 
[33] used this empirical relation to reduce the scatter in the Hubble diagram to 
a < 0.2 mag in V for a sample of nearly 30 SNIa from the CTSS search. 

Impressed by the success of the Z\mi5(B) parameter, Riess et al. [79] devel- 
oped the multi-color light curve shape method (MLCS), which parameterized the 
shape of SN light curves as a function of their absolute magnitude at maximum. 
This method also included a sophisticated error model and fitted observations 
in all colors simultaneously, allowing a color excess to be included. This color 
excess, which we attribute to intervening dust, enabled the extinction to be mea- 
sured. Another method that has been used widely in cosmological measurements 
with SNIa is the "stretch" method described in Perlmutter et al. [74,77]. This 
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method is based on the observation that the entire range of SNIa hght curves, 

at least in the B and V-bands, can be represented with a simple time stretching 
(or shrinking) of a canonical light curve. The coupled stretched B and V light 
curves serve as a parameterized set of light curve shapes [26], providing many 
of the benefits of the MLCS method but as a much simpler (and constrained) 
set. This method, as well as recent implementations of Z\mi5(B) [24,65], also 
allows extinction to be directly incorporated into the SNIa distance measure- 
ments. Other methods that correct for intrinsic luminosity differences or limit 
the input sample by various criteria have also been proposed to increase the 
precision of type la SNe as distance indicators [9,17,93,95]. While these latter 
techniques are not as developed as the Z\mi5(B), MLCS, and stretch methods, 
they all provide distances that are comparable in precision, roughly a = 0.18 
mag about the inverse square law, equating to a fundamental precision of SNIa 
distances of ^ 6% (0.12 mag), once photometric uncertainties and peculiar ve- 
locities are removed. Finally, a "poor man's" distance indicator, the snapshot 
method [80], combines information contained in one or more SN spectra with 
as little as one night's multi-color photometry. This method's accuracy depends 
critically on how much information is available. 



3 Cosmological Parameters 

The standard model for describing the global evolution of the Universe is based 
on two equations that make some simple, and hopefully valid, assumptions. If 
the Universe is isotropic and homogenous on large scales, the Robertson- Walker 
Metric, 



ds^ = dt^ - a{t) 



1 — kr^ 



(4) 



gives the line element distance(s) between two objects with coordinates r,0 and 
time separation, t. The Universe is assumed to have a simple topology such 
that, if it has negative, zero, or positive curvature, k takes the value —1,0,1, 
respectively. These models of the Universe are said to be open, flat, or closed, 
respectively. The dynamic evolution of the Universe needs to be input into the 
Robertson- Walker Metric by the specification of the scale factor a{t), which gives 
the radius of curvature of the Universe over time - or more simply, provides the 
relative size of a piece of space at any time. This description of the dynamics of 
the Universe is derived from General Relativity, and is known as the Friedman 
equation 



. id/af ='-^-^. (5) 



The expansion rate of our Universe (H), is called the Hubble parameter (or 
the Hubble constant, Hq, at the present epoch) and depends on the content of 
the Universe. Here we assume the Universe is composed of a set of components, 
each having a fraction, ]7j, of the critical density 
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with an equation of state which relates the density, p,, and pressure, p,, as Wi = 
Pi/ Pi- For example, Wi takes the value for normal matter, +1/3 for photons, 
and -1 for the cosmological constant. The equation of state parameter does not 
need to remain fixed; if scalar fields are present, the effective w will change over 
time. Most reasonable forms of matter or scalar fields have Wi > —1, although 
nothing seems manifestly forbidden. Combining Eqs. (4-6) yields solutions to 
the global evolution of the Universe [13]. 

The luminosity distance, Dl, which is defined as the apparent brightness 
of an object as a function of its redshift z - the amount an object's light has 
been stretched by the expansion of the Universe - can be derived from Eqs. (4- 
6) by solving for the surface area as a function of z, and taking into account 
the effects of time dilation [25,26,50,82] and energy dimimition as photons get 
stretched traveling through the expanding Universe. Dl is given by the numer- 
ically integrable equation, 



i5L = ^(l + ^)«o'^'^W rdz'[y2n,{l + z'f^^-^-^o{l + z'f]-'^'}. (7) 
^0 Jo Y 

^(a;) = sin(a;), x, or sinh(x) for closed, flat, and open models respectively, and 
the curvature parameter kq, is defined as kq = f2i — 1. 

Historically, Eq. (7) has not been easily integrated and has been expanded 
in a Taylor series to give 

where the deceleration parameter, qo, is given by 

90 = ^$^^»(l + 3wi). (9) 

i 

Prom Eq. (9) we can see that, in the nearby Universe, the luminosity distances 
scale linearly with redshift, with Hq serving as the constant of proportionality. 
In the more distant Universe, Dl depends to first order on the rate of accel- 
eration/deceleration (go) or, equivalently, on the amount and types of matter 
that make up the Universe. For example, since normal matter has wm = and 
the cosmological constant has wm = — 1, a universe composed of only these two 
forms of matter /energy has qo = i7M/2 — f^A- In a universe composed of these 
two types of matter, if i?yi < Om/'^, Qo is positive and the universe is decelerat- 
ing. These decelerating universes have Dl smaller as a function of z than their 
accelerating counterparts. 

If distance measurements are made at a low-^ and a small range of redshift 
at higher z (e.g., 0.3 > z > 0.5), there is a degeneracy between Hm and i?^. 
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Fig. 1. Dl expressed as distance modulus (m — M) for four relevant cosmological 
models; SIm = 0, f2A = (empty Universe, solid line); Om = 0.3, Qa ~ {short 
dashed line); Om — 0.3, Oa ~ 0.7 {hatched line); and Om ~ 1-0, Oa ~ {long dashed 
line). In the bottom panel the empty Universe has been subtracted from the other 
models to highlight the differences. 

It is impossible to pin down the absolute amount of either species of matter. 
One can only determine their relative dominance, which, at z = 0, is given by 
Eq. (9). However, Goobar and Perlmutter [27] pointed out that by observing 
objects over a larger range of high redshift (e.g., 0.3 > z > 1.0) this degeneracy 
can be broken, providing a measurement of the absolute fractions of f2M and 

To illustrate the effect of cosmological parameters on the luminosity distance, 
in Fig. 1 we plot a series of models for both A and non-il universes. In the top 
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Fig. 2. Dl for a variety of cosmological models containing Hm ~ 0.3 and Ox ~ 0.7 
with a constant (not time- varying) equation of state Wx- The Wx ~ model has been 
subtracted off to highlight the differences between the various models 



panel, the various models show the same linear behavior at z < 0.1 with models 
having the same Hq being indistinguishable to a few percent. By z = 0.5 the 
models with significant A are clearly separated, with luminosity distances that 
are significantly further than the zero-A universes. Unfortunately, two perfectly 
reasonable universes, given our knowledge of the local matter density of the 
Universe {Hm ^ 0.2), one with a large cosmological constant, Ha = 0.7, f2M = 
0.3 and one with no cosmological constant, f^M = 0.2, show differences of less 
than 25%, even to redshifts of z > 5. Interestingly, the maximum difference 
between the two models is at z 0.8, not at large z. Fig. 2 illustrates the effect 
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of changing the equation of state of the non-matter, dark energy component, 

assuming a flat universe, f2tot = 1. If wc arc to discern a dark energy component 
that is not a cosmological constant, measurements better than 5% are clearly 
required, especially since the differences in this diagram include the assumption 
of flatness and also fix the value of Hm- In fact, to discriminate among the full 
range of dark energy models with time varying equations of state will require 
much better accuracy than even this challenging goal. 

4 Meeisuring the Hubble Constant 

Schmidt et al. [86], using a sample of 16 SNII, estimated Hq = 73±6(statistical)±7 
(systematic) using EPM. This estimate is independent of other rungs in the cx- 
tragalactic distance ladder, the most important of which are the Cepheids, which 
currently calibrate most other distance methods (such as SNIa). The Cepheid 
and EPM distance scales, compared galaxy to galaxy, agree to within 5% and are 
consistent within the errors [16,52]. This provides confidence that both methods 
are providing accurate distances. 

The current nearby SNIa sample [24,32,41,84] contains more than 100 objects 
(Fig. 3), and accurately defines the slope in the Hubble diagram from < z < 0.1 
to 1%. To measure Hq, SNIa must still be externally calibrated with Cepheids, 
and this calibration is the major limitation to measuring Hq with SNIa. Two 
separate teams have analyzed the Cepheids and SNIa but have obtained di- 
vergent values for the Hubble constant. Saha et al. [88] find Hq = 59 ± 6, 
whereas Frecdman ct al. [20] find Hq ~ 71 ± 2 ± (6 systematic). Of the 12 
SNIa for which there are Cepheid distances to the host galaxy (SN1895B*, 
SN1937C*, SN1960F*, SN1972E, SN1974G*, SN1981B, SN1989B, SN1990N, 
SN1991T, SN1998cq, SN1998bu, and SN1999by), four were observed by non- 
digital means (marked by *) and are best excluded from analysis on the grounds 
that non-digital photometry routinely has systematic errors far greater than 0.1 
mag. Jha [41] has compared the SNIa distances using an updated version of 
MLCS to the Cepheid host galaxy distances measured by the two Hubble Space 
Telescope (HST) teams. Using only the digitally observed SNIa, he finds, using 
distances from the SNIa project of Saha et al. [88], Hq = 66±3± (7 systematic) 

km s~^ Mpc~^. Applying the same analysis to the Key Project distances by 
Preedman et al. [20] gives i?o = 76 ± 3 ± (8 systematic) km s~^ Mpc~^ This 
difference is not due to SNIa errors, but rather to the different ways the two 
teams have measured Cepheid distances with HST. The two values do overlap 
when the systematic uncertainties are included, but it is still uncomfortable that 
the discrepancies are so large, particularly when some systematic uncertainties 
are common between the two teams. 

At present, SNe provide the most convincing constraints with Hq ~ 70 ± 10 

km s""'^ Mpc""'^. However, future work on measuring Hq lies not with the SNe 
but with the Cepheid calibrators, or possibly in using other primary distance 
indicators such as EPM or the Sunyaev-Zeldovich effect. 
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Fig. 3. The Hubble diagram for SNIa from 0.01 > z > 0.2 [24,33,41,84]. The 102 
objects in this range have a residual about the inverse square hne of ~ 10%. 



5 The Measurement of Acceleration 

The intrinsic brightness of SNIa allow them to be discovered to z > 1.5 with 
current instrumentation (while a comparably deep search for type II SNe would 
only reach redshifts of z 0.5). In the 1980s, however, finding, identifying, and 
studying even the impressively luminous type la SNe was a daunting challenge, 
even towards the lower end of the redshift range shown in Fig. 1. At these 
redshifts, beyond z ~ 0.25, Fig. 1 shows that relevant cosmological models could 
be distinguished by differences of order 0.2 mag in their predicted luminosity 
distances. For SNIa with a dispersion of 0.2 mag, 10 well observed objects should 
provide a 3cr separation between the various cosmological models. It should be 
noted that the uncertainty described above in measuring Hq is not important 
in measuring the parameters for different cosmological models. Only the relative 
brightness of objects near and far is being exploited in Eq. (7) and the absolute 
value of Hq scales out. 

The first distant SN search was started by the Danish team of N0rgaard- 
Nielsen et al. [57]. With significant effort and large amounts of telescope time 
spread over more than two years, they discovered a single SNIa in a z = 0.3 
cluster of galaxies (and one SNII at z = 0.2) [35,57]. The SNIa was discovered 
well after maximum light on an observing night that could not have been pre- 
dicted, and was only marginally useful for cosmology. However, it showed that 
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such high redshift SNe did exist and could be found, but that they would be 

very difficult to use as cosmological tools. 

Just before this first discovery in 1988, a search for high redshift type la 
SNe using a then novel wide field camera on a much larger (4m) telescope was 
begun at the Lawrence Berkeley National Laboratory (LBNL) and the Center 
for Particle Astrophysics, at Berkeley. This search, now known as the Supernova 
Cosmological Project (SCP), was inspired by the impressive studies of the late 
1980s indicating that extremely similar type la SN events could be recognized 
by their spectra and light curves, and by the success of the LBNL fully robotic 
low-redshift SN search in finding 20 SNe with automatic image analysis [56,67]. 

The SCP targeted a much higher redshift range, z > 0.3, in order to measure 
the (presumed) deceleration of the Universe, so it faced a different challenge 
than the CTSS search. The high redshift SNe required discovery, spectroscopic 
confirmation, and photometric follow up on much larger telescopes. This precious 
telescope time could neither be borrowed from other visiting observers and staff 
nor applied for in sufficient quantities spread throughout the year to cover all 
SNe discovered in a given search field, and with observations early enough to 
establish their peak brightness. Moreover, since the observing time to confirm 
high redshift SNe was significant on the largest telescopes, there was a clear 
"chicken and egg" problem: telescope time assignment committees would not 
award follow-up time for a SN discovery that might, or might not, happen on a 
given run (and might, or might not, be well past maximum) and, without the 
follow-up time, it was impossible to demonstrate that high redshift SNe were 
being discovered by the SCP. 

By 1994, the SCP had solved this problem, first by providing convincing ev- 
idence that SNc, such as SN1992bi, could be discovered near maximum (and 
K-corrected) out to z = 0.45 [73], and then by developing and successfully 
demonstrating a new observing strategy that could effectively guarantee SN dis- 
coveries on a predetermined date, all before or near maximum light [70-72,76]. 
Instead of discovering a single SN at a time on average (with some runs not 
finding one at all), the new approach aimed to discover an entire "batch" of 
half-a-dozen or more type la SNc at a time by observing a much larger mimber 
of galaxies in a single two or three day period a few nights before new Moon. By 
comparing these observations with the same observations taken towards the end 
of dark time almost three weeks earlier, it was possible; to select just those SNe 
that were still on the rise or near maximum. The chicken and egg problem was 
solved, and now the follow-up spectroscopy and photometry could be applied 
for and scheduled on a pre-specified set of nights. The new strategy worked ~ 
the SCP discovered batches of high redshift SNe, and no one would ever again 
have to hunt for high-redshift SNe without the crucial follow-up scheduled in 
advance. 

The High-Z SN Search (HZSNS) was conceived at the end of 1994, when this 
group of astronomers became convinced that it was both possible to discover 
SNIa in large numbers at z > 0.3 by the efforts of Perlmutter et al. [70-72], and 
also use them as precision distance indicators as demonstrated by the CTSS 
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group [32]. Since 1995, the SCP and HZSNS have both worked avidly to obtain 
a significant set of high redshift SNIa. 

5.1 Discovering SNIa 

The two high redshift teams both used this pre-scheduled discovery and follow- 
up batch strategy. They each aimed to use the observing resources they had 
available to best scientific advantage, choosing, for example, somewhat different 
exposure times or filters. 

Quantitatively, type la SNe are rare events on an astronomer's time scale 
- they occur in a galaxy like the Milky Way a few times per millennium (see, 
e.g., [12,60,61] and the chapter by Cappellaro in this volume). With modern 
instruments on 4 meter-class telescopes, which observe 1/3 of a square degree to 
R = 24 mag in less than 10 minutes, it is possible to search a million galaxies to 
z < 0.5 for SNIa in a single night. 

Since SNIa take approximately 20 days to rise from undetectable to max- 
imum light [81], the three- week separation between observing periods (which 
equates to 14 rest frame days at z = 0.5) is a good filter to catch the SNe on 
the rise. The SNe are not always easily identified as new stars on the bright 
background of their host galaxies, so a relatively sophisticated process must be 
used to identify them. The process, which involves 20 Gigabytes of imaging data 
per night, consists of aligning a previous epoch, matching the image star profiles 
(through convolution), and scaling the two epochs to make the two images as 
identical as possible. The difference between these two images is then searched 
for new objects which stand out against the static sources that have been largely 
removed in the differencing process [73,74,76,87]. The dramatic increase in com- 
puting power in the 1980s was an important element in the development of this 
search technique, as was the construction of wide-field cameras with ever larger 
CCD detectors or mosaics of such detectors [104]. 

This technique is very efficient at producing large numbers of objects that 
are, on average, at or near maximum light, and does not require unrealistic 
amounts of large telescope time. It does, however, place the burden of work on 
follow-up observations, usually with different instruments on different telescopes. 
With the large number of objects discovered (50 in two nights being typical), 
a new strategy is being adopted by both the SCP and HZSNS teams, as well 
as additional teams like the Canada France Hawaii Telescope (CFHT) legacy 
survey, where the same fields are repeatedly scanned several times per month, in 
multiple colors, for several consecutive months. This type of observing program 
provides both discovery of new objects and their follow up, all integrated into one 
efficient program. It does require a large block of time on a single telescope - a 
requirement which was not politically feasible in years past, but is now possible. 

5.2 Obstacles to Measuring Luminosity Distances at High-Z 

As shown above, the distances measured to SNIa are well characterized a.t z < 
0.1, but comparing these objects to their more distant counterparts requires great 
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care. Selection effects can introduce systematic errors as a function of redshift, 

as can imccrtain K-corrcctions and a possible evolution of the SNIa progenitor 
population as a function of look-back time. These effects, if they are large and 
not constrained or corrected, will limit our ability to accurately measure relative 
luminosity distances, and have the potential to reduce the efficacy of high-^; type 
la SNe for measuring cosmology [74,77,83,87]. 

K-Corrections: As SNe are observed at larger and larger redshifts, their light 

is shifted to longer wavelengths. Since astronomical observations are normally 
made in fixed band passes on Earth, corrections need to be applied to account for 
the differences caused by the spectrum shifting within these band passes. These 
corrections take the form of integrating the spectrum of an SN over the relevant 
band passes, shifting the SN spectrum to the correct redshift, and re-integrating. 
Kim et al. [43] showed that these effects can be minimized if one does not use 
a single bandpass, but instead chooses the bandpass closest to the rcdshiftcd 
rest- frame bandpass, as they had done for SN1992bi [73]. They showed that the 
inter-band K-correction is given by 



Kij{z) = 2.5log 



i^ + z] 



J F{\)S.,{X)d\ JZ{\)Sji\)d\ 
J F{X/{1 + z))Sj{X)d\) J Z{X)S,{\)dX 



(10) 



where Kij{z) is the correction to go from filter i to filter j, and Z{X) is the 
spectrum corresponding to zero magnitude of the filters. 

The brightness of an object expressed in magnitudes, as a function of z is 



mi{z) = 5 log 



Mpc 



■25 + M^+K,j{z), 



(11) 



where Dl{z) is given by Eq. (7), Mj is the absolute magnitude of object in 
filter j, and Kij is given by Eq. (10). For example, for Hq = 70 km s~^ Mpc~^, 
and Dl = 2835 Mpc [Qm = 0.3, /?yi = 0.7), at maximum light a SNIa has 
Mb = —19.5 mag and a Kbb, = —0.7 mag. We therefore expect an SNIa at 
z = 0.5 to peak at mR ~ 22.1 mag for this set of cosmological parameters. 
K-correction errors depend critically on three uncertainties: 

1. Accuracy of spectrophotometry of SNe. To calculate the K-correction, the 
spectra of SNe are integrated in Eq. (10). These integrals are insensitive 
to a grey shift in the flux calibration of the spectra, but any wavelength 
dependent flux calibration error will translate into erroneous K-corrections. 

2. Accuracy of the absolute calibration of the fundamental astronomical stan- 
dard systems. Eq. (10) shows that the K-corrections are sensitive to the 
shape of the astronomical band passes and to the zero points of these band 
passes. 

3. Accuracy of the choice of SNIa spectrophotometry template used to calculate 
the corrections. Although a relatively homogenous class, there are variations 
in the spectra of SNIa. If a particular object has, for example, a stronger 
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calcium triplet than the average SNIa, the K-corrections will be in error 
unless an appropriate subset of SNIa spectra are used in the calculations. 

The first error should not be an issue if correct observational procedures are 
used on an instrument that has no fundamental problems. The second error is 
currently estimated to be small (~ 0.01 mag), based on the consistency of spec- 
trophotometry and broadband photometry of the fundamental standards, Sirius 
and Vega [5]. To improve this uncertainty will require new, careful experiments 
to accurately calibrate a star, such as Vega or Sirius (or a White Dwarf or solar 
analog star), and to carefully infer the standard bandpass that defines the pho- 
tometric system in use at telescopes. The third error requires a large database to 
match as closely as possible an SN with the spectrophotometry used to calculate 
the K-corrections. Nugent et al. [58] have shown that extinction and color are 
related and, by correcting the spectra to force them to match the photometry 
of the SN needing K-corrections, that it is possible to largely eliminate errors 1 
and 3, even when using spectra that are not exact matches (in epoch or in fine 
detail) to the SNIa being K-corrected. Scatter in the measured K-corrections 
from a variety of telescopes and objects allows us to estimate the combined size 
of the eff'(;ct for the first and third errors. These appear to be ^ 0.01 mag for 
redshifts where the high- 2; and low- 2; filters have a large region of overlap (e.g., 
R-band matched to B-band at z = 0.5). 

Extinction: In the nearby Universe we see SNIa in a variety of environments, 
and about 10% have significant extinction [30]. Since we can correct for extinc- 
tion by observing two or more wavelengths, it is possible to remove any first 
order eflfects causcid by a changing average extinction of SNIa as a function of z. 
However, second order effects, such as possible evolution of the average proper- 
ties of intervening dust, could still introduce systematic errors. This problem can 
also be addressed by observing distant SNIa over a decade or so of wavelength 
in order to measure the extinction law to individual objects. Unfortunately, this 
is observationally very expensive. Current observations limit the total system- 
atic effect to < 0.06 mag, as most of our current data is based on two color 
observations. 

An additional problem is the existence of a thin veil of dust around the Milky 
Way. Measurements from the Cosmic Background Explorer (COBE) satellite 
accurately determined the relative amount of dust around the Galaxy [89], but 
there is an uncertainty in the absolute amount of extinction of about 2 — 3%. 
This uncertainty is not normally a problem, since it affects everything in the sky 
more or less equally. However, as we observe SNe at higher and higher redshifts, 
the light from the objects is shifted to the red and is less affected by the Galactic 
dust. Our present knowledge indicates that a systematic error as large as 0.06 
mag is attributable to this uncertainty. 

Selection Effects: As we discover SNe, we are subject to a variety of selection 
effects, both in our nearby and distant searches. The most significant effect is the 
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Malmquist Bias - a selection effect which leads magnitude limited searches to 

find brighter than average objects near their distance hmit since brighter objects 
can be seen in a larger volume than their fainter counterparts. Malmquist Bias 
errors are proportional to the square of the intrinsic dispersion of the distance 
method, and because SNIa arc such accurate distance indicators these errors are 
quite small, ~ 0.04 mag. Monte Carlo simulations can be used to estimate such 
selection effects, and to remove them from our data sets [74,76,77,87]. The total 
uncertainty from selection effects is ^ 0.01 mag and, interestingly, may be worse 
for lower redshift objects because they are, at present, more poorly quantified. 

Gravitational Lensing: Several authors have pointed out that the radiation 

from any object, as it traverses the large scale structure between where it was 
emitted and where it is detected, will be weakly lensed as it encounters fluctu- 
ations in the gravitational potential [37,45,100]. On average, most of the light 
travel paths go through under-dense regions and objects appear de-magnified. 
Occasionally, the light path encounters dense regions and the object becomes 
magnified. The distribution of observed fluxes for sources is skewed by this pro- 
cess such that the vast majority of objects appear slightly fainter than the canon- 
ical luminosity distance, with the few highly magnified events making the mean 
of all light paths unbiased. Unfortunately, since we do not observe enough ob- 
jects to capture the entire distribution, unless we know and include the skewed 
shape of the lensing a bias will occur. At z = 0.5, this lensing is not a significant 
problem: If the Universe is flat in normal matter, the large scale structure can 
induce a shift of the mode of the distribution by only a few percent. However, 
the effect scales roughly as z^, and by z = 1.5 the effect can be as large as 25% 
[38] . While corrections can be derived by measuring the distortion of background 
galaxies near the line of sight to each SN, at z > 1, this problem may be one 
which ultimately limits the accuracy of luminosity distance measurements, un- 
less a large enough sample of SNe at each redshift can be used to characterize 
the lensing distribution and average out the effect. For the z ^ 0.5 sample, the 
error is < 0.02 mag, but it is much more significant at z > 1 (e.g., for SN1997ff) 
[4,55], especially if the sample size is small. 

Evolution: SNIa are seen to evolve in the nearby Universe. Hamuy et al. [29] 
plotted the shape of the SN light curves against the type of host galaxy. SNe 
in early hosts (galaxies without recent star formation) consistently show light 
curves which rise and fade more quickly than SNe in late-type hosts (galaxies 
with on-going star formation). However, once corrected for light curve shape the 
luminosity shows no bias as a function of host type. This empirical investigation 
provides reassurance for using SNIa as distance indicators over a variety of stel- 
lar population ages. It is possible, of course, to devise scenarios where some of 
the more distant SNe do not have nearby analogues, so as supernovae are stud- 
ied at increasingly higher redshifts it can become important to obtain detailed 
spectroscopic and photometric observations of every distant SN to recognize and 
reject examples that have no nearby analogues. 
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In principle, it should be possible to use differences in the spectra and light 

curves between nearby and distant SNe, combined with theoretical modeling, to 
correct any differences in absolute magnitude. Unfortunately, theoretical inves- 
tigations are not yet advanced enough to precisely quantify the effect of these 
differences on the absolute magnitude. A different, empirical approach to handle 
SN evolution [10] is to divide the SNe into subsamples of very closely matched 
events, based on the details of the their light curves, spectral time series, host 
galaxy properties, etc. A separate Hubble diagram can then be constructed for 
each subsample of SNe, and each will yield an independent measurement of the 
cosmological parameters. The agreement (or disagreement) between the results 
from the separate subsamples is an indicator of the total effect of evolution. A 
simple, first attempt at this kind of test has been performed by comparing the 
results for SNe found in elliptical host galaxies to SNe found in late spirals or 
irregular hosts, and the cosmological results from these subsamples were found 
to agree well [91]. 

Finally, it is possible to move to higher redshifts and see if the SNe deviate 

from the predictions of Eq. (7). At a gross level, we expect an accelerating Uni- 
verse to be decelerating in the past because the matter density of the Universe 
increases with redshift, whereas the density of any dark energy leading to ac- 
celeration will increase at a slower rate than this (or not at all in the case of 
a cosmological constant). If the observed acceleration is caused by some sort of 
systematic effect, it is likely to continue to increase (or at least remain steady) 
with z, rather than disappear like the effects of dark energy. A first comparison 
has been made with SN1997ff at z ~ 1.7 [85], and it seems consistent with a 
decelerating Universe at that epoch. More objects are necessary for a definitive 
answer, and these should be provided by several large programs that have been 
discovering such type la SNe at the W.M. Keck Telescope I (KECK 1), Subaru 
Telescope), and ffST telescopes. 

5.3 High Redshift SNIa Observations 

The SCP [74] in 1997 presented their first results with 7 objects at a redshift 
around z = 0.4. These objects hinted at a decelerating Universe with a measure- 
ment of J?M = O.SS^Q gQ, but were not definitive. Soon after, the SCP published 
a further resuff, with a z ~ 0.84 SNIa observed with the KECK I and HST 
added to the sample [75] , and the HZSNS presented the results from their first 
four objects [22,87]. The results from both teams now ruled out a /?m = 1 Uni- 
verse with greater than 95% significance. These findings were again superceded 
dramatically when both teams announced results including more SNe (10 more 
HZSNS SNe, and 34 more SCP SNe) that showed not only were the SN obser- 
vations incompatible with a f2M = 1 Universe, they were also incompatible with 
a Universe containing only normal matter [77,83]. Fig. 4 shows the Hubble dia- 
gram for both teams. Both samples show that SNe are, on average, fainter than 
would be expected, even for an empty Universe, indicating that the Universe is 
accelerating. The agreement between the experimental results of the two teams 



Measuring cosmology with Supernovae 17 




r . — . — — ] 

^ 0.01 0.10 1.00 

z 

Fig. 4. Upper panel: The Hubble diagram for high redshift SNIa from both the HZSNS 
[83] and the SCP [77]. Lower panel: The residual of the distances relative to a Hm — 0.3, 
Oa = 0.7 Universe. The z < 0.15 objects for both teams are drawn from CTSS sample 
[32], so many of these objects are in common between the analyses of the two teams. 
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Fig. 5. The confidence regions for both HZSNS [83] and SCP [77] for Qm, Oa- The two 
experiments show, witli remarkable consistency, that Qa > is required to reconcile 
observations and theory. The SCP result is based on measurements of 42 distant SNIa. 
(The analysis shown here is uncorrected for host galaxy extinction;see [77] for the 
alternative analyses with host extinction correction, which is shown to make little 
difference in this data set.) The HZSNS result is based on measurements of 16 SNIa, 
including 6 snapshot distances [80], of which two arc SCP SNe from the 42 SN sample. 
The z < 0.15 objects used to constrain the fit for both teams are drawn from the CTSS 
sample [32], so many of these objects are common between the analyses by the two 
teams. 
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Fig. 6. Left panel; Contours of Qm versus Wx from current observational data. Right 
Panel: Contours of Qm versus Wx from current observational data, where the current 
value of Qm is obtained from the 2dF redshift survey. For both panels Qm + f?!, = 1 
is taken as a prior. 



is spectacular, especially considering the two programs have worked in almost 
complete isolation from each other. 

The easiest solution to explain the observed acceleration is to include an ad- 
ditional component of matter with an equation of state parameter more negative 
than w < —1/3; the most familiar being the cosmological constant {w — — 1). 
Fig. 5 shows the joint confidence contours for values of and f^A from both 
experiments. If we assume the Universe is composed only of normal matter and 
a cosmological constant, then with greater than 99.9% confidence the Universe 
has a non-zero cosmological constant or some other form of dark energy. 

Of course, we do not know the form of dark energy which is leading to the 
acceleration, and it is worthwhile investigating what other forms of energy are 
possible additional components. Fig. 6 shows the joint confidence contours for 
the HZSNS+SCP observations for and Wx (the equation of state of the 
unknown component causing the acceleration) . Because this introduces an extra 
parameter, we apply the additional constraint that fiM + ^x = 1, as indicated by 
the CMB experiments [14] . The cosmological constant is preferred, but anything 
with a. w < —0.5 is acceptable [23,77]. Additionally, we can add information 
about the value of Qm , as supplied by recent 2dF redshift survey results [98] , as 
shown in the 2nd panel, where the constraint strengthens to w < —0.6 at 95% 
confidence [69]. 




Fig. 7. Future expected constraints on dark energy: Left panel: Estimated 68% con- 
fidence regions for a constant equation of state parameter for tlie dark energy, w, 
versus mass density, for a ground-based study witli 200 SNe between z = 0.3 — 0.7 
(open contours), and for the satellite-based SNAP experiment with 2,000 SNe between 
z = 0.3 — 1.7 (filled contours). Both experiments are assumed to also use 300 SNe 
between z = 0.02 — 0.08. A flat cosmology is assumed (based on Cosmic Microwave 
Background (CMB) constraints) and the inner (solid line) contours for each experiment 
include tight constraints (from large scale structure surveys) on i?M, at the ±0.03 level. 
For the SNAP experiment, systematic uncertainty is taken as dm = 0.02(z/1.7), and 
for the ground-based experiment, dm = 0.03(^/0.5). Such ground-based studies will 
test the hypothesis that the dark energy is in the form of a cosmological constant, 
for which w = —1 at all times. Middle panel: The same confidence regions for the 
same experiments not assuming the equation of state parameter, w, to be constant, 
but instead marginalizing over w', where w{z) = wo + w'z. (Weller and Albrecht [103] 
recommend this parameterization of w{z) over the others that have been proposed to 
characterize well the current range of dark energy models.) Note that these planned 
ground-based studies will yield impressive constraints on the value of w today, wo , even 
without assuming constant w. In fact, these constraints are comparable to the current 
measurements of w assuming it is constant (shown in the right panel of Fig. 6). Right 
panel: Estimated 68% confidence regions of the first derivative of the equation of state, 
w' , versus its value today, wo, for the same experiments. 



6 The Future 

How far can we push the SN measurements? Finding more and more SNe allows 
us to beat down statistical errors to arbitrarily small levels but, ultimately, 
systematic effects will limit the precision to which SNIa magnitudes can be 
applied to measure distances. Our best estimate is that it will be possible to 
control systematic eflFects from ground-based experiments to a level of ~ 0.03 
mag. Carefully controlled ground-based experiments on 200 SNc will reach this 
statistical uncertainty in z = 0.1 redshift bins in the range z = 0.3 — 0.7, and 
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is achievable within five years. A comparable quality low redshift sample, with 
300 SNe in z = 0.02 — 0.08, will also be achievable in that time frame [2]. 

The SuperNova/ Acceleration Probe (SNAP) collaboration^ has proposed to 
launch a dedicated cosmology satellite [1,68] - the ultimate SNIa experiment. 
This satellite will, if funded, scan many square degrees of sky, discovering well 
over a thousand SNIa per year and obtain their spectra and light curves out to 
z = 1.7. Besides the large numbers of objects and their extended redshift range, 
space-based observations will also provide the opportunity to control many sys- 
tematic effects better than from the ground [21,53]. Fig. 7 shows the expected 
precision in the SNAP and ground-based experiments for measuring iv, assuming 
a flat Universe. Perhaps the most important advance will be the first studies of 
the time variation of the equation of state w (see the right panel of Fig. 7 and 
[40,103]). 

With rapidly improving CMB data from interferometers, the satellites Mi- 
crowave Anisotropy Probe (MAP) and Planck, and balloon-based instrumenta- 
tion planned for the next several years, CMB measurements promise dramatic 
improvements in precision on many of the cosmological parameters. However, the 
CMB measurements are relatively insensitive to the dark energy and the epoch 
of cosmic acceleration. SNIa are currently the only way to directly study this 
acceleration epoch with sufficient precision (and control on systematic uncer- 
tainties) that we can investigate the properties of the dark energy, and any time 
dependence in these properties. This ambitious goal will require complementary 
and supporting measurements of, for example, Qm from CMB, weak lensing, 
and large scale structure. The SN measurements will also provide a test of the 
cosmological results independent from these other techniques, which have their 
own systematic errors. Moving forward simultaneously on these experimental 
fronts offers the plausible and exciting possibility of achieving a comprehensive 
measurement of the fundamental properties of our Universe. 
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